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1  Introduction 


Classical  military  terrain  analysis  typically  includes  an  interpretation  and 
evaluation  of  numerous  geophysical  and  biophysical  characteristics  of  the  earth’s 
surface  within  some  specific  area  of  interest.  Geophysical  parameters  include 
surface  materials  (including  both  soils  and  surface  geology),  surface  configura¬ 
tion  (i.e.,  slope  and  aspect),  surface  drainage  (i.e.,  landform  class),  and  water  re¬ 
sources.  The  primary  biophysical  parameter  is  a  detailed  description  of  the 
unique  vegetation  cover  types.  The  interaction  and  spatial  dependencies  of  adja¬ 
cent  or  overlapping  geophysical  and  biophysical  parameters  can  result  in  highly 
complex  terrain  that  potentially  requires  extensive  reconnaissance  to  accurately 
and  precisely  define  and  delineate.  Remotely  sensed  data  sourqes  may  offer  the 
most  cost-effective  method  for  large  area  terrain  analysis.  However,  most  air¬ 
borne  and  space  imagery  acquired  for  both  manual  and  automated  (or  semi- 
automated)  feature  extraction  is  acquired  from  overhead  orientation.  Therefore, 
with  the  sensor  oriented  nadir  to  the  earth’s  surface,  identifying  and  delineating 
physical  phenomena  that  are  subordinate  to  (i.e.,  beneath)  taller  and  wider 
dominant  feature  classes  may  be  highly  inaccurate,  if  not  unattainable. 

Terrain  analysis  data  base  specifications  list  desired  “overlay”  information  for  a 
variety  of  geophysical  parameters  (DMA,  1982).  Several  critical  vegetation 
structural  and  compositional  variables  are  listed  as  well,  including: 

1.  Vegetation  type 

2.  Canopy  closure  percentage 

3.  Stem  spacing 

4.  Vegetation  height 

5.  Vegetation  roughness  factors 

6.  Tree  stem  diameter 

7.  Undergrowth. 

Each  of  these  vegetation  characteristics  is  identified  as  an  essential  element  of 
terrain  information  (HQDA,  1990)  to  be  delineated  onto  hardcopy  map  overlays. 
Each  overlay  is  the  product  of  some  level  of  detailed  manual  interpretation  of  one 
or  more  remotely  sensed  images.  The  monoscopic  and  stereoscopic  surface  fea¬ 
ture  classification  relies  on  the  interpreter’s  ability  to  accurately  identify  the 
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dominant  landscape  theme  (e.g.,  forest  type,  shrub  community,  grass  species  as¬ 
sociation)  and  then  to  accurately  infer  the  subordinate  landscape  themes  (e.g., 
understory  vegetation,  soil  conditions,  water  class). 

This  research  effort  examines  the  empirical  relationship  between  the  species 
composition  and  physical  structure  of  the  dominant  vegetation  layer  and  the  bio¬ 
physical  density  of  the  subordinate  (i.e.,  understory)  vegetation  layers.  Current 
military  terrain  analysis  methods  employ  subjective  procedures  for  estimating 
the  nature  and  extent  of  understory  vegetation  strata.  Tactical  data  base  specifi¬ 
cations  used  by  the  National  Intelligence  and  Mapping  Agency  (NIMA)  assign 
understory  vegetation  density  estimates  to  only  two  possible  classes:  (1)  greater 
than  50  percent  understory  density,  and  (2)  less  than  50  percent  density  or  unde¬ 
termined.  These  extremely  broad  class  definitions  provide  only  limited  input  to 
the  various  tactical  terrain  models,  including:  cross  country  mobility,  cover  and 
concealment,  line-of-sight,  and  bivouac  areas.  Furthermore,  the  characteristic  of 
“50  percent  understory  density”  can  be  interpreted  in  a  variety  of  ways  and 
therefore  does  not  represent  a  truly  empirical  measure  of  the  biomass  of  woody 
and  herbaceous  vegetation  below  the  dominant  woody  overstory.  Therefore,  there 
is  a  need  to: 

•  improve  on  the  definitions  of  understory  and  understory  density  relevant  to 
tactical  military  terrain  data  base  generation  protocols 

•  enhance  the  analyst’s  ability  to  accurately  and  precisely  estimate  understory 
vegetation  physical  and  compositional  characteristics  through  the  interpreta¬ 
tion  of  the  dominant  overstory  characteristics 

•  allow  the  data  to  be  quickly  and  accurately  predicted  using  areal  photogra¬ 
phy  and  high-resolution  multi-spectral  imagery. 
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2  Objectives 

The  ability  to  observe,  record,  and  quantify  dominant  vegetation  features  over 
large  areas  has  been  developed  and  refined  over  the  past  60  years  following  the 
widespread  use  of  both  large-  and  small-scale  aerial  photography  and  imagery. 
In  the  last  30  years,  digital  technologies  have  been  used  to  create  vegetation  type 
maps.  However,  the  ability  to  directly  observe  subordinate  vegetation  strata  re¬ 
mains  a  predictive  procedure  for  typical  passive,  optical  remote  sensing  devices. 
Some  applications  with  active  sensors,  such  as  airborne  radar  and  lidar  systems, 
have  attempted  to  directly  quantify  understory  vegetation  densities  with  mar¬ 
ginal  accuracy. 

Therefore,  the  purpose  of  this  project  is  to  begin  to  develop  some  preliminary 
predictive  relationships  between  forest  overstory  characteristics  and  forest  un¬ 
derstory  densities.  The  specific  objectives  include: 

•  to  develop  a  field  sampling  strategy  that  ensures  accurate  and  precise  quan¬ 
tification  of  both  forest  overstory  and  woody  understory  species  composition 
and  vertical  structure 

•  to  develop  predictive  relationships  using  forest  overstory  parameters  as  the 
independent  variables  and  woody  understory  stem  density  as  the  dependent 
variable. 
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3  Background 

The  study  of  the  relative  quantities  of  vegetation  biomass  within  multi-layered 
or  multi-storied  forest  types  has  been  applied  to  a  wide  variety  of  ecological  ar¬ 
eas.  Examples  of  studies  that  have  investigated  the  differences  in  the  horizontal 
distributions  of  the  vertically  arranged  composition  and  structure  of  multi¬ 
layered  forest  canopies,  include: 

•  wildlife  habitat  evaluations 

•  forest  successional  research 

•  vertical  species  diversity  studies  (especially  in  tropical  and  neotropical  re¬ 
gions). 

The  characteristics  that  define  each  individual  layer  in  a  multi-tiered  forest  can¬ 
opy  differ  drastically  from  one  forest  type  to  another.  Past  site  disturbances, 
whether  natural  (e.g.,  fire,  flood,  storm  event)  or  anthropogenic  (land  clearing, 
controlled  bum,  pollution  impacts)  likely  play  the  greatest  role  in  determining 
the  current  forest  stand  dynamic  conditions.  After  the  impacts  from  disturbance, 
local  and  regional  landscape  characteristics  typically  control  the  number  of  spe¬ 
cies  within  a  specific  forest  type  and  their  vertical  placement  in  the  perennial 
vegetated  canopy.  Some  of  the  site-specific  factors  associated  with  landform  that 
directly  influence  forest  understory  development  include:  climate,  landform  posi¬ 
tion  (e.g.,  slope,  aspect  and  exposure),  and  soil  moisture. 

This  study  collected  forest  stand  dynamics  data  within  mid-latitude,  temperate, 
mixed  (deciduous  and  conifer)  ecosystems.  The  sample  sites  have  been  impacted 
by  manmade  disturbances  for  roughly  350  years.  Two  of  the  study  sites,  how¬ 
ever,  show  little  anthropogenic  impacts  over  the  last  250  to  300  years.  As  com¬ 
pared  to  some  other  mid-latitude  temperate  forest  regions  in  other  parts  of  the 
world,  particularly  central  and  southern  Europe,  this  period  of  human  distur¬ 
bance  is  very  brief. 

The  typical  forests  sampled  in  the  study  exhibited  either  mixed  hardwood  or 
mixed  hardwood/conifer  species  composition.  Stand  structures  were  typically 
representative  of  uneven-aged  canopies  with  multi-layered  canopies  that  in¬ 
cluded  distinct  dominant  stems  and  one  or  more  layers  of  intermediate  and  sup¬ 
pressed  stems.  The  understory  layers  ranged  from  a  few  widely  scattered  stems 
to  fairly  dense  stands  of  seedlings  and  saplings.  The  exceptions  to  these  repre- 
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sentative  hardwood  and  mixed  stands  were  pine  plantations.  Where  managed, 
the  understory  component  at  these  sites  was  typically  minimal.  When  left  un¬ 
managed,  the  understory  density  of  these  plantations  was  substantial. 
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4  Methodology 


Field  Data  Collection  Sites 

Sample  sites  were  selected  based  on  prior  forest  sampling  knowledge  and  experi¬ 
ence  within  and  around  Central  Virginia.  While  an  attempt  was  made  to  ran¬ 
domly  select  plot  locations  consideration  was  given  to  ease  of  access.  In  addition, 
an  effort  was  made  to  obtain  samples  from  a  wide  variety  of  forest  types  (young 
and  old-age  stands,  variable  species  composition,  etc.).  Nineteen  field  data  col¬ 
lects  were  taken  from  14  October  1999  through  July  2000  in  Central  Virginia. 
Data  were  collected  at  56  plot  sites  spread  out  over  seven  locations  (Figure  1). 

The  field  data  collection  sites  were  located  in  Northern  and  North  Central  Vir¬ 
ginia  on  the  eastern  fringes  of  the  Piedmont.  The  topography  of  the  areas  gener¬ 
ally  consists  of  flat  fields  to  gently  rolling  slopes.  Local  relief  is  generally  less 
than  100  ft  (30.5  m).  In  a  number  of  data  collect  areas,  the  landscape  was  more 
dissected  as  streams  cut  through  forested  areas.  Soils  of  the  region  are  generally 
Ultisols  and  Vertisols. 

The  plot  locations  and  respective  number  of  plots  at  each  location  were: 

1.  FortAP.  Hill,  Bowling  Green,  VA.  Fort  A.P.  Hill  is  a  U.S.  Army  training  installa¬ 
tion  located  approximately  15  to  20  mi  (1  mi  =  1.61  km)  south  of  Fredericksburg, 
VA.  The  topography  is  generally  flat,  with  numerous  bottomlands  and  marshy 
areas.  There  are  also  locations  on  the  Fort  that  present  a  more  dissected  land¬ 
scape  where  creeks  cut  through  forested  areas  (37  plots). 

2.  Topographic  Engineering  Center  (TEC),  Alexandria,  VA.  TEC  is  situated  ap¬ 
proximately  20  mi  SSW  of  Washington,  DC  in  southern  Fairfax  County  The  to¬ 
pography  surrounding  TEC  is  generally  flat  with  numerous  marshy  and  swampy 
areas,  many  of  them  seasonal  in  nature  (three  plots). 

3.  Caledon  Natural  Area,  King  George,  VA.  This  location  is  a  special  area  set  aside 
by  the  State  of  Virginia  and  managed  by  Virginia’s  Department  of  Conservation 
and  Recreation.  It  is  located  approximately  25  mi  east  of  Fredericksburg,  VA. 

The  land,  a  donation  of  a  wealthy  family,  consists  of  forested  hills  and  wetlands 
adjacent  to  the  Potomac  River.  The  forests  on  the  property  are  described  as  sec¬ 
ond-growth  stands  (four  plots). 


4.  Montpelier,  VA.  Montpelier  is  located  approximately  40  mi  west  southwest  of 
Fredericksburg,  VA.  It  is  the  former  estate  of  President  James  Madison.  (The 
grounds  are  managed  by  an  historical  foundation.)  The  topography  is  generally 
rolling.  Numerous  areas  on  the  estate  contain  individual  old-growth  specimens 
and  small,  undisturbed  old-growth  stands  (four  plots). 

5.  Fredericksburg  Battlefield,  Fredericksburg  VA.  This  battlefield  is  located  ap¬ 
proximately  1  mi  from  downtown  Fredericksburg,  VA  and  is  maintained  by  the 
National  Park  Service.  The  topography  is  flat  to  gently  rolling.  Pre-Civil  War 
photographs  of  this  site  show  an  abundance  of  agricultural  fields  with  minimal 


Figure  1 .  Data  Collection  Sites. 


©1 995  by  Ray  Sterner,  Applied  Physics  Laboratory,  The  Johns  Hopkins  University.  Map  used  with  permission 
Available  at:  fermi.jhuapl.edu/states. 


woody  vegetation  (two  plots). 


6.  Spotsylvania  National  Military  Park,  Fredericksburg  VA.  This  site  is  located  ap¬ 
proximately  10  mi  southwest  of  Fredericksburg,  VA.  The  topography  is  flat  to 
gently  rolling.  As  with  the  Fredericksburg  Battlefield,  this  area  was  composed  of 
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agricultural  fields  during  the  Civil  War.  What  woody  vegetation  that  existed  was 
cut  down  by  troops  for  fuel  and  fortifications  (two  plots) 

7.  Chancellorsville  Battlefield,  Fredericksburg,  VA.  This  site  is  located  approxi¬ 
mately  10  mi  west  of  Fredericksburg,  VA.  The  topography  is  primarily  flat.  Its 
histoxy  mirrors  that  of  the  previous  two  sites  (four  plots). 


Plot  Design 

Fixed-area,  circular  plots  were  designed  to  accurately  and  precisely  sample  both 
the  overstory  and  understory  layers.  Experience  from  previous  forest  sampling 
projects  suggested  a  plot  radius  of  10  to  15  m.  To  minimize  sampling  efforts 
within  each  plot  and  to  maximize  the  total  number  of  plots  a  10-m  radius  (~33  ft) 
was  selected.  Wenger  (1984)  recommends  that,  for  forests  in  the  eastern  United 
States,  sample  plots  should  range  from  1/50  to  1/25  hectare  (ha).  The  10-m  ra¬ 
dius  plot  is  equivalent  to  a  plot  size  of  approximately  1/32-ha.  Furthermore, 
given  that  this  work  has  an  ultimate  application  in  forest  attribute  extraction 
from  remotely  sensed  imagery,  a  20-m  diameter  plot  theoretically  adequately 
characterizes  both  a  SPOT  4  and  Landsat  ETM  7  pixel. 


Data  Collection  Procedures 
Overstory  Stems 

A  plot  center  was  established  at  a  representative  location  within  the  selected 
forest  stand.  A  small  wire  flag  was  placed  in  the  ground  to  mark  the  center  of 
the  sampling  area.  Plot  coordinates  were  obtained  using  a  PLGR  (GPS  receiver). 
Flagging  was  then  positioned  at  10  m  from  the  plot  center  at  each  of  the  4  cardi¬ 
nal  compass  directions.  Working  clockwise  from  a  dominant  tree  within  the  plot, 
the  following  attributes  were  recorded  for  each  dominant  and  codominant  over¬ 
story  tree  using  standard  measurement  techniques  (see  Avery  and  Burkhart, 
1994): 

1.  Distance  from  plot  center  (meters) 

2.  Azimuth  (true  north)  using  plot  center  as  the  reference  (degrees) 

3.  Species 

4.  Diameter  at  breast  height  (cm) 

5.  Total  height 

6.  Canopy  closure 
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7.  Crown  diameter 

8.  Crown  canopy  class. 

Overstory  individuals  were  defined  as  all  woody  stems  with  a  diameter  at  breast 
height  (d.b.h.)  greater  than  4  in.  (~10  cm).  The  distance  and  azimuth  measures 
were  collected  to  quantitatively  calculate  average  minimum  distance  between 
the  stems.  This  measure  represents  an  estimate  of  stem  spacing  within  the  plot. 
Tree  species  was  recorded.  The  d.b.h.  measurements  were  collected  using  either 
a  caliper  or  diameter  tape.  Trunk  diameter  is  consistently  measured  at  ap¬ 
proximately  4.5  ft  (1.37  m)  above  the  ground.  Precision  was  within  0.1  in  (~0.25 
cm).  Tree  height  was  measured  using  a  clinometer.  Precision  was  on  the  order  of 
1  ft  (~0.3  m). 

Canopy  class  describes  the  vertical  position  of  each  overstory  crown  with  respect 
to  the  surrounding  individuals.  For  statistical  analysis  in  this  study,  the  over¬ 
story  was  deemed  to  be  only  the  dominant  and  codominant  stems.  The  crowns  of 
these  stems  can  be  readily  seen  by  airborne  platforms.  Intermediate  and  sup¬ 
pressed  stems,  not  viewable  from  above,  were  merged  with  the  understory  stems 
during  analysis. 

Definitions  of  each  of  the  four  crown  canopy  classes  is  provided  by  Smith  (1962): 

1.  Dominant :  Trees  with  crowns  extending  above  the  general  level  of  the  crown 
cover  and  receiving  frill  light  from  above  and  partially  from  the  side,  larger  than 
the  average  trees  in  the  stand,  and  with  crowns  well  developed,  but  possibly 
somewhat  crowded  on  the  sides. 

2.  Codominant :  Trees  with  crowns  forming  the  general  level  of  the  crown  cover  and 
receiving  frail  light  from  above,  but  comparatively  little  from  the  sides,  usually 
with  medium-sized  crowns  more  or  less  crowded  on  the  sides. 

3.  Intermediate :  Trees  shorter  than  the  dominant  and  codominant  with  crowns 
pushing  into  the  crown  cover,  receiving  little  light  from  above  and  none  from  the 
sides,  usually  with  small  crowns  considerably  crowded  on  the  sides. 

4.  Suppressed :  Also  known  as  overtopped,  these  trees' have  crowns  that  fall  entirely 
below  the  general  level  of  the  crown  cover,  receiving  no  direct  light  either  from 
above  or  from  the  sides 

For  all  stems  identified  as  either  dominant  or  codominant,  crown  diameter  was 
measured  and  recorded.  Crown  diameter  was  not  measured  for  crowns  classified 
as  intermediate  or  suppressed.  Two  measurements  were  obtained,  one  repre¬ 
senting  the  major  crown  axis  and  the  other  the  minor  crown  axis.  Crown  axis  is 
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determined  by  estimating  the  drip-line  of  the  crown.  These  values  were  then 
summed  and  divided  by  2  to  obtain  an  average  crown  diameter.  Precision  was  on 
the  order  of  1  ft  (0.3  m). 

Ancillary  plot  attributes  recorded  at  each  plot  location  included: 

1.  Slope  (percent) 

2.  Aspect  (degrees) 

3.  Landform  type 

4.  Canopy  closure  (percent) 

5.  General  soil  description 

6.  Description  of  overall  shrub/sapling,  seedling,  and  herbaceous  cover. 

Understory  Stems 

For  the  purpose  of  this  study,  understory  stems  were  defined  as  all  woody  indi¬ 
viduals  with  d.b.h.  less  than  4  in  (-10  cm)  and  total  height  greater  than  6  ft  (~1.8 
m).  Using  these  criteria,  measured  understory  specimens  consisted  of  saplings 
and  taller  shrubs.  Species,  height,  crown  diameter  and  d.b.h.  were  measured  for 
all  sampled  understory  stems.  (Preliminary  data  analysis  suggested  that  under- 
stoiy  height  and  crown  diameter  showed  practically  no  correlation  with  any 
other  stand  variables.  Thus,  these  two  measurements  were  not  recorded  on 
roughly  one-third  of  the  sample  sites.) 

Plant  material  that  was  less  than  6  ft  (-1.8  m)  in  total  height  consisted  of  tree 
seedlings,  dwarf  shrubs,  and  the  herbaceous  layer  (grasses,  forbs,  and  ferns).  No 
direct  measurements  of  this  vegetation  stratum  were  recorded.  However,  a  de¬ 
scription  of  these  nonsampled  understory  vegetation  members  was  noted  in  the 
plot  logs. 


Preliminary  Data  Management 

All  of  the  acquired  plot  information  was  entered  onto  an  Excel  spreadsheet  and 
imported  into  Statistica  (StatSoft,  1999)  for  manipulation  and  analysis.  Sum¬ 
mary  statistics  were  then  generated  for  each  plot,  including  means,  standard  de¬ 
viations,  and  variances  for  overstory  and  understory  height,  crown  diameter,  and 
d.b.h. 
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Overstory  and  understory  stem  densities  were  derived  and  reported  as  the  total 
number  of  estimated  stems  per  hectare  (to  convert  to  stems  per  acre,  divide  by 
2.471).  Mean  minimum  distance  between  the  dominant  and  codominant  stems 
was  reported  for  each  plot.  Finally,  a  description  of  the  species  composition  of 
each  plot  was  used  to  classify  the  sample  sites  into  broad  vegetation  types,  in¬ 
cluding: 

1.  Pure  hardwood 

2.  Mixed  hardwood/pine 

3.  Mixed  pine/hardwood 

4.  Pure  pine. 


Statistical  Analysis 

Due  to  the  paucity  of  previous  work  on  this  subject  found  in  the  literature,  the 
authors  had  no  a  priori  hypotheses  concerning  the  relationships  between  the 
overstory  characteristics  and  understory  stem  density.  Therefore,  all  measured 
variables  (stem  characteristics  and  plot  attributes)  were  included  in  the  analysis. 
The  initial  approach  to  uncovering  relationships  was  through  employment  of  ex¬ 
ploratory  data  analysis  (EDA)  techniques.  EDA  is  an  approach  to  data  analysis 
that  postpones  the  usual  assumptions  about  what  kind  of  model  the  data  follow 
with  the  more  direct  approach  of  allowing  the  data  itself  to  reveal  its  underlying 
structure  and  model.  Most  EDA  techniques  are  graphical  in  nature.  Graphics 
provide  analysts  with  open-minded  exploration  using  their  intuitive  pattern- 
recognition  capabilities.  Primary  EDA  techniques  include  scatterplots,  histo¬ 
grams,  residual  plots,  probability  plots,  and  plots  of  simple  statistics  such  as 
means  and  standard  deviations.  The  data,  therefore,  are  used  to  suggest  the  ap¬ 
propriate  model(s)  that  fit  the  data  itself. 
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5  Results 


Summary  Statistics 

A  total  of  4314  individual  trees  were  measured — 1051  overstory  stems  and  3263 
understory  stems — across  56  individual  plots.  Table  1  lists  the  total  number  of 
overstory  stems  sampled  by  species.  Loblolly  pine  ( Pinus  taeda )  was  the  most 
frequently  encountered  softwood  species  (22.4  percent),  while  white  oak  ( Quercus 
alba )  was  the  most  frequently  sampled  hardwood  species  (18.8  percent).  A  total 
of  22  unique  deciduous  species  were  encountered  in  the  overstory  while  only  five 
different  conifer  species  were  found  across  all  56  plots.  Table  2  lists  the  number 
of  understory  stems  tallied  by  species.  American  holly  (Ilex  opaca ),  one  of  three 
evergreen  hardwood  species  commonly  found  in  North  America,  was  the  most 
frequently  sampled  understory  species  (17.5  percent).  Highbush  blueberry  ( Vac - 
cinium  corymbosum )  and  dogwood  ( Cornus  florida )  were  the  second  and  third, 
respectively,  most  frequently  encountered  understory  species.  Virginia  pine 
( Pinus  virginiana),  loblolly  pine  (P.  taeda),  and  eastern  red  cedar  ( Juniperus  vir- 
giniana )  were  the  only  conifer  species  occupying  the  under  story.  Their  small 
numbers  can  be  attributed  to  the  fact  that  they  are  all  highly  shade  intolerant. 

Table  3  displays  summary  univariate  overstory  statistics  for  each  plot.  Plots  9d 
and  12c,  both  pine  plantations,  support  the  greatest  density  of  overstory  stems  at 
923  and  924  per  hectare,  respectively.  The  plots  with  the  fewest  number  of  stems 
per  hectare  included  19b  (64  stems  per  hectare),  plot  5d  (95  stems  per  hectare), 
and  plot  7a  (96  stems  per  hectare).  Each  of  these  sites  was  oak  dominated. 

Mean  stem  diameter  (d.b.h.)  was  calculated  for  each  plot.  The  sites  with  the 
greatest  average  stem  diameters  included  plots  19a-d.  These  oak  dominated 
stands  averaged  in  diameter  between  18.2  and  32.2  in.  Plots  12c  and  9a  dis¬ 
played  the  small  average  stem  diameters  with  5.1  to  6.4  in.,  respectively. 

The  greatest  stand  densities  as  measured  by  basal  area  are  found  in  plots  19a- 
19d,  with  estimated  total  basal  areas  per  hectare  of  12.9  to  16.3  sq  ft  (1.2  to  1.5 
m2).  These  sites  all  support  oak  dominated  overstories.  Also,  plots  4a  and  4b, 
older  second-growth  stands  in  Caledon  Natural  Area,  have  basal  areas  of  14.7 
and  16.6  sq  ft  (1.37  and  1.54  m2)  respectively. 
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Table  1.  Overstory  Stems — Species  and  Counts. 


|  Common  Name 

Latin  Name 

Stem  Count 

Pinus  taeda 

235 

22.4 

white  oak 

Quercus  alba 

198 

18.8 

Virginia  pine 

Pinus  virginiana 

105 

10.0 

yellow  poplar 

Liriodendron  tulipifera 

104 

9.9 

Liquidambar  styraciflua 

70 

6.7 

|  southern  red  oak 

Quercus  falcata 

55 

5.2 

Acer  rubrum 

50 

4.8 

beech 

Fagus  grandifolia 

43 

4.1 

hickory 

Carya  tomentosa 

40 

3.8 

black  oak 

Quercus  velutina 

34 

3.2 

Nyssa  sylvatica 

23 

2.2 

|  chestnut  oak 

Quercus  prinus 

20 

1.9 

Populus  grandidentata 

18 

1.7 

post  oak 

scarlet  oak 

northern  red  oak 

dogwood 
pitch  pine 

river  birch 

persimmon 

sassafras 

longleaf  pine 

sycamore 

cherry 

willow  oak 

blackjack  oak 

eastern  red  cedar 

Quercus  stellata 

Quercus  coccinea 

Quercus  rubra 

Comus  florida 

Pinus  rigida 

Betula  nigra 

Diospyros  virginiana 

Sassafras  albidum 

Pinus  palustris 

Plantanus  occidentalis 

Prunus  serotina 

Quercus  phellos 

Quercus  marilandica 

Juniperus  virginiana 

56 

5.3 

Totals 

1051 

100.0 

The  sites  supporting  the  lowest  basal  areas  included  4c,  with  only  1.97  sq  ft  (0.18 
m2),  18a,  with  3.01  sq  ft  (0.28  m2)  and  12c,  with  4.26  sq  ft  (0.4  m2).  The  first  two 
plots,  4c  and  18a,  were  hardwood  dominated,  while  12c  was  a  young,  unmanaged 
pine  plantation. 

Plots  19  (a-d)  and  4  (a-c)  maintained  the  overall  tallest  stands,  with  mean  overstory 
heights  ranging  from  103  to  124  ft  (31.4  to  37.8  m).  The  shortest  stands  sampled 
were  in  plot  12c,  a  young,  unmanaged  pine  plantation  with  an  average  total  height 
of  only  27.3  ft  (8.3  m),  and  plot  9a,  another  pine  dominated  site  supporting  stems 
only  41.3  ft  (12.6  m)  tall  on  average.  Plots  16a-c  and  17a-d  were  pine  plantations 
with  average  overstoiy  canopy  heights  from  54.0  to  77.9  ft  (23.7  m). 
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Table  2.  Understory  Steins — Species  and  Counts. 


Common  Name 

Latin  Name 

Stem  Count 

ESBS9HI 

holly 

Ilex  opaca 

570 

17.5 

Vaccinium  corymbosum 

380 

11.6 

Comus  florida 

376 

11.5 

I  southern  red  oak 

Quercus  falcata 

365 

11.2 

Liquidambar  styraciflua 

337 

10.3 

beech 

Fagus  grandifolia 

210 

6.4 

white  oak 

Quercus  alba 

200 

6.1 

hickory 

Carya  tomentosa 

194 

5.9 

Acer  rubrum 

136 

4.2 

1  black  oak 

Quercus  velutina 

70 

2.1 

Nyssa  sylvatica 

64 

2.0 

Populus  grandidentata 

57 

1.7 

mountain  laurel 

Kalmia  latifolia 

53 

1.6 

sassafras 

Sassafras  albidum 

50 

1.5 

Liriodendron  tulipifera 

32 

1.0 

Pinus  virginiana 

32 

1.0 

magnolia 

Magnolia  Virginia 

25 

0.8 

chestnut  oak 

Quercus  prinus 

20 

0.6 

cherry 

Pmnus  serotina 

18 

0.6 

loblolly 

Pinus  taeda 

18 

H9H 

hornbeam 

Carpinus  caroliniana 

14 

0.4 

willow  oak 

Quercus  phellos 

11 

0.3 

eastern  red  cedar 

Juniperus  virginiana 

9 

0.3 

spicebush 
downy  serviceberry 
devil's  walking  stick 
striped  maple 
American  chestnut 

Lindera  benzoin 

Amelanchier  arborea 
Aralia  spinosa 

Acer  pensylvanicum 
Castanea  dentata 

22 

0.7 

Totals 

3263 

100.0 

The  final  overstory  parameter  of  importance  is  average  crown  diameter.  The 
largest  average  crown  diameters  were  measured  in  4b  and  4c,  both  with  61.3  ft 
(18.7  m)  on  average.  These  plots  were  dominated  by  sweetgum  and  oak,  respec¬ 
tively,  and  were  located  at  Caledon  Natural  Area,  an  old,  second-growth  forest. 
The  smallest  crown  diameters  were  found  in  plot  9a  at  11.5  ft  (3.5  m).  Many 
more  plots  possessed  an  average  crown  diameter  between  12  and  15  ft  (3.7  and 
4.6  m)  wide.  All  the  small  crown  plots  were  established  within  pine  plantations. 


Table  4  presents  summary  univariate  understory  statistics  by  plot.  As  stated 
above,  the  number  of  understory  variables  collected  in  each  plot  was  reduced  a 
little  more  than  half  way  through  the  sampling. 


Table  3.  Overstory  Summary  Statistics. 
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Therefore,  the  only  parameters  that  were  consistently  sampled  within  all  56  sites 
were  the  density  of  understory  stems  and  d.b.h.  The  minimum  estimated  density 
of  568  stems  per  hectare  was  recorded  in  plot  4b  which  maintained  a  sweet  gum 
overstory.  The  understory  species  in  plot  4b  were  dominated  by  American  holly. 
The  greatest  understory  density  was  found  in  plot  12c  with  an  estimate  6210 
understory  stems  per  hectare  under  a  pure  pine  canopy  in  a  young  unmanaged 
plantation.  The  dominant  understory  species  in  this  plantation  site  was  south¬ 
ern  red  oak.  As  a  whole,  these  summary  statistics  suggest  an  adequate  repre¬ 
sentation  to  perform  a  preliminary  investigation  of  the  potential  relationships 
between  forest  understory  density  to  overstory  parameters. 


Exploratory  Data  Analyses 

EDA  techniques  revealed  a  number  of  interesting  relationships  and  structures 
between  overstory  characteristics  and  understory  stem  density.  More  classical 
statistical  techniques  were  then  employed.  Classification  techniques  were  used 
to  try  to  organize  the  data  into  meaningful  groupings.  Hierarchical  cluster 
analysis  (joining-tree  clustering)  and  nonhierarchical  clustering  (k-means)  were 
employed  to  uncover  natural  clusters  or  groups  in  the  data.  Discriminant  analy¬ 
sis  was  then  used  in  an  effort  to  determine  the  most  important  variables  that 
discriminate  between  groups.  These  techniques,  it  was  hoped,  would  divide  the 
groups  into  those  with  like  overstory  characteristics  enabling  the  examination  of 
the  corresponding  understory  stem  densities  in  each  group.  The  aforementioned 
techniques  were  also  performed  on  subsets  of  the  data.  The  56  plots  were  di¬ 
vided  by  observable  break-points  in  the  height,  crown,  and  d.b.h.  distributions. 
The  plots  were  additionally  divided  by  species  and  combinations  of  species.  This 
proved  less  than  successful  primarily  because  of  the  small  number  of  data  points 
used  in  the  subset  analyses. 

The  results  of  these  analyses  determined  that  using  summarized  data  for  all  56 
plots  was  the  best  approach  in  developing  preliminary  predictive  relationships. 
Then  both  simple  and  multiple  regression  were  performed  on  the  plot  data. 


Regression  Results 

The  following  discussion  presents  the  results  of  both  simple  and  multiple  regres¬ 
sion  in  estimating  understory  stem  density.  For  the  purposes  of  this  study,  un¬ 
derstory  stems  were  limited  to  those  with  d.b.h.  values  >1.0  in  (>2.54  cm)  and 
height  >6.0  ft  (~1.8  m).  They  also  included  the  intermediate  and  suppressed 
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overstory  stems,  which  cannot  be  viewed  from  above.  This  excludes  many  of  the 
multi-stemmed  shrubs.  In  all  cases,  simple  regression  took  the  form  of  a  single 
term  function  with  intercept.  Examined  first  were  the  relationships  between  the 
measured  parameters  of  the  overstory  itself.  Next,  all  single  overstory  parame¬ 
ters  are  regressed  against  understory  stem  density.  Finally,  multiple  regression 
was  used  to  estimate  the  understory  stem  density. 

Overstory  Structural  Relationships — Single  Predictor  (Independent) 

Variable 

The  first  focus  was  on  the  overstory  itself.  This  is  what  can  be  directly  observed 
and  measured  by  the  analyst.  As  shown  by  the  scatterplot  patterns  in  Figures  2 
through  9,  the  overstory  exhibits  a  high  degree  of  structure.  If  the  overstory 
structure  were  purely  random,  then  the  relationships  between  the  measured 
overstory  and  the  understory  might  very  well  be  nonexistent.  Regressions  were 
performed  using  Version  4.0  of  TableCurve  2D  (SPSS,  1997).  This  software  proc¬ 
esses  x,y  data  points  with  over  8000  linear  and  nonlinear  equations  to  determine 
a  best-fit.  It  was  decided  to  limit  the  body  of  potential  equations  to  single  and 
two  term  linear  functions.  Hence,  each  of  the  best-fit  equations  appearing  below 
was  selected  from  a  body  of  over  1600  linear  equations. 

In  addition  to  the  form  of  the  best-fit  equation  and  the  “a”  and  “b”  coefficients, 
each  graph  contains  a  number  of  summary  statistics.  The  coefficient  of  determi¬ 
nation,  r2,  indicates  the  proportion  of  variability  in  the  dependent  variable  that  is 
explained  by  the  independent  (predictor)  variable.  The  Adj  (adjusted)  r2  at¬ 
tempts  to  correct  the  r2  to  more  closely  reflect  the  goodness  of  the  fit  of  the  model. 
It  incorporates  the  number  of  cases  and  number  of  independent  variables.  The 
standard  error  of  the  estimate  (S.E.E.)  is  a  measure  of  the  dispersion  of  the  ob¬ 
served  values  about  the  regression  line.  The  “F”  value  (Fstat)  shows  how  well 
the  regression  model  fits  the  data.  If  the  probability  associated  with  this  value  is 
small,  one  can  reject  the  hypothesis  that  the  r2  value  =  0. 

Figure  2  shows  the  relationship  between  overstory  stem  density  and  overstory 
mean  height  for  the  56  forest  plots.  Older  forests  are  those  that  appear  in  the 
upper  left-hand  portion  of  the  graph.  Younger  forests  are  shorter  and  have 
higher  overstory  stem  densities  and  appear  in  the  lower  right-hand  portion  of 
the  graph.  The  relationship  exhibits  a  moderate  r2  value  of  ~0.69. 

Figure  3  shows  the  relationship  between  the  overstory  stem  density  per  hectare 
and  overstory  mean  crown  diameter  for  the  56  forest  plots.  Older  forests  appear 
in  the  upper  left-hand  comer  of  the  graph  and  younger  forests  in  the  lower  right- 
hand  comer  of  the  graph.  The  relationship  exhibits  a  moderate  r2  of  0.72. 


Mean  Crown  Diameter  (ft.)  I  “  Overstory  Mean  Height  (ft.) 
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Overstory  Stem  Density  vs.  Overstory  Mean  Height 
y=a+bx°-5 

r2=0. 68789867  Adj  r2=0.67612126  S.E.E  =12.872786  Fstat=1 19.02073 

a=1 43. 15934 
b=-3.1 766972 


re  2.  Overstory  Stem  Density  vs.  Overstory  Mean  Height. 


Overstory  Stem  Density  vs.  Overstory  Mean  Crown  Diameter 

y=a+fc»£5 

r2=0.72214633  Adj  r2=0.71 166129  S.E.E.=6.3399841  Fstat=140.3469 
a=57. 352672 
b-1. 6989537 


Figure  3.  Overstory  Stem  Density  vs.  Overstory  Mean  Crown  Diameter. 
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Figure  4  shows  the  relationship  between  the  overstory  stem  density  per  hectare 
and  overstory  mean  d.b.h.  for  the  56  forest  plots.  Older  forests  appear  in  the  up¬ 
per  left-hand  corner  of  the  graph  and  younger  forests  in  the  lower  right-hand 
comer  of  the  graph.  The  relationship  exhibits  a  moderate  r2  of  ~0.77. 

Figure  5  shows  the  relationship  between  the  overstory  stem  density  per  hectare 
and  the  mean  minimum  distance  of  the  dominant  and  codominant  stems  for  the 
56  forest  plots.  Older  forests  appear  in  the  upper  left-hand  comer  of  the  graph 
and  younger  forests  in  the  lower  right-hand  comer  of  the  graph. 

Figure  6  shows  the  relationship  between  the  overstory  mean  height  and  the 
overstory  mean  crown  for  the  56  forest  plots.  Older  forests  appear  in  the  upper 
right-hand  comer  of  the  graph  and  younger  forests  in  the  lower  left-hand  comer 
of  the  graph. 

Figure  7  shows  the  relationship  between  the  overstory  mean  height  and  the 
overstory  mean  d.b.h.  for  the  56  forest  plots.  Older  forests  appear  in  the  upper 
right-hand  comer  of  the  graph  and  younger  forests  in  the  lower  left-hand  comer 
of  the  graph. 

Figure  8  shows  the  relationship  between  the  overstory  mean  crown  diameter  and 
the  overstory  mean  d.b.h.  for  the  56  forest  plots.  Older  forests  appear  in  the  up¬ 
per  right-hand  comer  of  the  graph  and  younger  forests  in  the  lower  left-hand 
comer  of  the  graph.  The  r2  is  0.87,  the  highest  for  all  of  the  overstory  relation¬ 
ships. 

The  results  shown  in  Figure  8  mirror  those  of  Krause,  Puffenberger,  and  Camp¬ 
bell  (1999)  in  their  study  of  crown  diameter-d.b.h.  relationships.  It  is  interesting 
to  note  that  the  1999  study  only  considered  trees  with  round,  uniform  crowns, 
whereas  this  study  incorporated  all  overstory  trees  regardless  of  crown  uniform¬ 
ity.  The  simple  linear  relationship  still  holds  true. 

Figure  9  shows  the  relationship  between  the  overstory  mean  height  and  the 
mean  minimum  distance  of  the  dominant  and  codominant  stems  for  the  56  forest 
plots.  Older  forests  appear  in  the  upper  right-hand  comer  of  the  graph  and 
younger  forests  in  the  lower  left-hand  comer  of  the  graph.  This  relationship  is 
marginal  as  the  residuals  are  highly  heteroskedastic  (fanning  out  when  the 
mean  height  >  60  ft  [18.3  m])  and  the  r2  is  quite  poor  (0.40). 
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Overstory  Stem  Density  vs.  Overstory  Mean  d.b.h. 

y0-5=3+bx0.5 

r2=0.77529716  Adj  r2=0.76665475  S.E.E.=2.3103283  Fstat=1 82.86707 


a=5. 5006076 
b=-0.101 38018 


4.  Overstory  Stem  Density  vs.  Overstory  d.b.h. 


Overstory  Stem  Density  vs.  Mean  Minimum  Distance 
y=a+blnx 

r2=0.69982088  Adj  r2=0.68849337  S.E.E.=3.2455159  Fstat=1 25.89259 


a=52. 34461 
b=-6. 9864498 


Figure  5.  Overstory  Stem  Density  vs.  Mean  Minimum  Distance. 
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Overstory  Mean  Height  vs.  Overstory  Mean  Crown  Diameter 

y=a+bx2 

r2=0.69142801  Adj  r2=0.67978378  S.E.E.=6.6812597  Fstat=1 20.99968 


a=6. 9954 508 
b=0.0025396849 
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6.  Overstory  Mean  Height  vs.  Overstory  Mean  Crown  Diameter. 

Overstory  Mean  Height  vs.  Overstory  Mean  d.b.h. 
y=a+bx2 

r2=0. 82265577  DF  Adj  r2=0.81 596354  FitStdErr=2.0818032  Fstat=250.49257 


a=4. 8162727 
b=0.001 1385876 
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Figure  7.  Overstory  Mean  Height  vs.  Overstory  Mean  d.b.h. 


Minimum  Distance  (ft.)  ‘g  Overstory  Mean  d.b.h.  (in.) 


32 


ERDC/TEC  TR-01  -1 


Overstory  Mean  Crown  Diameter  vs.  Overstory  mean  d.b.h. 

y=a+bx 

r2=0.87936199  Adj  r2=0.87480961  S.E.E.=1.7788309  Fstat=393.62011 


a=3. 165458 
b=0. 39929725 
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re  8.  Overstory  Mean  Crown  Diameter  vs.  Overstory  Mean  d.b.h. 


Overstory  Mean  Height  vs.  Mean  Minimum  Distance 
y=a+bx 

r2=0.4006812  Adj  r2=0.3780654  S.E.E.=4.5858775  Fstat=36. 102297 


a=-0. 90261 054 
b=0. 16273023 


Figure  9.  Overstory  Mean  Height  vs.  Mean  Minimum  Distance. 
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Overstory/Understory  Relationships— Single  Predictor  (Independent) 

Variable 

All  of  the  measured  and  derived  overstory  and  plot  variables  were  used  indi¬ 
vidually  as  independent  variables  for  the  regression.  Table  5  contains  the  results 
of  the  regressions  using  single  independent  (predictor)  variable  to  estimate  un¬ 
derstory  stem  density  (stems  per  hectare).  The  r2  values,  even  in  the  best  cases, 
are  relatively  low  (~0.50).  Approximately  one-half  of  the  variability  in  the  un¬ 
derstory  stem  density  can  be  explained  by  the  overstory  variables.  Overstory 
mean  crown  diameter  is  the  best  predictor  of  understory  stem  density.  The  stan¬ 
dard  error  of  the  estimate  (value  =  520)  indicates  that  95  percent  of  all  under¬ 
story  stem  density  estimates  should  fall  within  ±  1040  understory  stems.  This 
translates  into  ±1  stem  for  every  103.5  sq  ft  (9.6  m2).  Figures  10  through  14 
show  plots  of  the  first  5  of  these  regressions  based  on  the  coefficient  of  determi¬ 
nation  (r2)  and  accompanying  Standard  Error  of  the  Estimate  (S.E.E.). 

Figure  10  shows  the  capability  of  overstory  mean  crown  diameter  in  estimating 
understory  stem  density.  Older  forests  appear  in  the  lower  right-hand  comer  of 
the  graph  and  younger  forests  in  the  upper  left-hand  corner  of  the  graph. 

Figure  11  shows  the  over  story  stem  density  per  hectare  as  a  predictor  of  under¬ 
story  stem  density.  Older  forests  appear  in  the  lower  left  portion  of  the  graph 
and  younger  stands  in  the  upper  right  comer 

Figure  12  shows  the  overstory  mean  height  as  a  predictor  of  understory  stem 
density.  Younger  stands  appear  in  the  upper  left-hand  portion  of  the  graph  and 
older  stands  appear  in  the  lower  right  hand  comer. 


Table  5.  Regression  Results  of  Understory  Stem  Density — Single  Predictor  Variable. 


Predictor  Variable 

1 

Ad¬ 

justed 

r2 

Standard 

Error  of 

the  Esti¬ 
mate 

Best-Fit 

Equation 

Regression  Coefficients 

As  predictor 

variable  increases 

the  number  of 

understory  stems 

a 

b 

OS  Mean  Crown  Diameter  (ft) 

0.56 

0.54 

520 

y  =  a+b/x2 

802.8 

292273.0 

Decreases 

OS  Stem  Density  -ha 

0.54 

0.52 

534 

y  =  a  +  bx2 

1023.4 

0.0026 

Increases 

OS  Mean  Height  (ft) 

0.51 

0.49 

547 

y  =  a  +  bx0.5 

-1832.6 

30027 

Decreases 

OS  Mean  D.B.H.  (in) 

0.47 

0.47 

560 

y  =  a  +  b/x2 

839.5 

81533.6 

Decreases 

Mean  Min.  Dist.  (ft)  of  OS  Stems 

0.42 

0.40 

596 

y  =  a  +  b/x 

291.3 

12779.7 

Decreases 

Variance  of  OS  Crown 

0.37 

0.34 

624 

y  =  a  +  b/x0.5 

727.0 

4288.0 

Decreases 

Variance  of  OS  Height 

0.20 

0.17 

702 

y  =  a  +  b/x2 

1300.0 

82683.0 

Decreases 

Slope  of  Plot  (%) 

0.22 

0.19 

730 

y  =  a  +  be-x 

1215.4 

988.0 

Decreases 

Variance  of  OS  D.B.H. 

0.13 

0.11 

732 

y  =  a  +  b/x0.5 

943.2 

1537.0 

Decreases 

Basal  Area  of  the  OS  (in2) 

0.02 

o 

o 

776 

y  =  a  +  bx3 

1535.8 

2.5E-08 

Decreases 
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Overstory  Mean  Crown  Diameter  vs.  Understory  Stem  Density 

y=a+b h£ 

r2=0. 561 52638  Adj  r2=0.54398744  S.E.E  =519.89958  Fstat=65.312585 


a=802. 83748 
b=292273.41 


Figure  10.  Overstory  Mean  Crown  Diameter  vs.  Understory  Stem  Density. 


Overstory  Stem  Density  vs.  Understory  Stem  Density 

y=a+bx2 

r2=0.5358253  Adj  r2=0.51725831  S.E.E.=534.91954  Fstat=58.872424 


a=1 023.4277 
b=0. 0025697131 


Figure  11.  Overstory  Stem  Density  vs.  Understory  Stem  Density. 
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Overstory  Mean  Height  vs.  Understory  Stem  Density 
y=a+b/>£-5 

r2=0. 51 287823  DF  Adj  r2=0.49339336  FitStdErr=547.98226  Fstat=53.696614 

a=-1 832.5727 
b=30027.23 
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Figure  12.  Overstory  Mean  Height  vs.  Understory  Stem  Density. 


Figure  13  shows  the  overstory  mean  d.b.h.  as  a  predictor  of  understory  stem  den¬ 
sity.  Younger  stands  appear  in  the  upper  left  comer  of  the  graph  and  older 
stands  in  the  lower  right  portion  of  the  graph. 

Figure  14  shows  the  mean  minimum  distance  of  the  overstory  stems  as  a  predic¬ 
tor  of  understory  stem  density.  Younger,  more  dense  stands  appear  in  the  upper 
left  of  the  graph  and  older,  less  dense  stands  appear  in  the  lower  right. 

Overstory/Understory  Relationships — Multiple  Regression  Result 

Multiple  regression  was  then  performed  on  the  data  set.  The  general  purpose  of 
multiple  regression  is  to  analyze  the  relationship  between  several  predictor  (in¬ 
dependent)  variables  (overstory  measures)  and  a  dependent  variable  (understory 
stem  density  per  hectare).  Three  different  types  of  multiple  regression  were  per¬ 
formed.  Standard  multiple  regression  enters  all  the  predictor  variables  at  once. 
Forward  stepwise  multiple  regression  adds  or  deletes  individual  independent 
variables  until  the  best  regression  model  is  obtained.  Backward  stepwise  multi¬ 
ple  regression  begins  with  all  independent  variables  in  the  model  and  then  re¬ 
moves  one  variable  at  a  time  until  the  best  model  is  obtained. 


Understory  Stems  per  hectare  c  Understory  Stems  per  hectare 
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Overstory  Mean  d.b.h.  vs.  Understory  Stem  Density 
y=a+b/>£ 

r2=0.48963895  Adj  r2=0.46922451  S.E.E.=560.90136  Fstat=48.929257 

a=839.45461 
b=81 533.608 


-e  13.  Overstory  Mean  d.b.h.  vs.  Understory  Stem  Density. 


Mean  Minimum  Distance  of  Overstory  Stems  vs.  Understory  Stem  Density 

y=a+b/x 

r2=0.42449147  Adj  r2=0.401471 13  S.E.E.=595.62595  Fstat=37.61728 
a=293. 11191 
b=1 2779.709 


Figure  14.  Mean  Minimum  Distance  of  Overstory  Stems  vs.  Understory  Stem  Density. 


Understory  Stems  per  hectare  Understory  Stems  per  hectare 
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Each  of  these  different  approaches  in  multiple  regression  was  performed  with 
both  a  y-intercept  present  and  with  the  y-intercept  set  to  zero.  For  variables  to 
be  accepted  in  the  equation,  their  p-values  had  to  be  less  than  0.05.  This  value, 
of  course,  is  arbitrary,  but  small  p-values  help  ensure  that  only  statistically  sig¬ 
nificant,  nonzero  coefficients  are  used  in  the  model. 

Figure  15  shows  the  results  for  the  multiple  regression  for  the  56  plots.  Ob¬ 
served  values  of  understory  stem  density  appear  on  the  y-axis  and  predicted  val¬ 
ues  along  the  x-axis.  The  best-fit  equation  had  no  y-intercept  and  was: 

Understory  stem  density  per  hectare  =  (2.997  x  over  stems  per  hectare)  + 

(27.788  x  overstory  mean  d.b.h.) 

The  multiple  regression  did  not  perform  as  well  as  when  single  predictor  vari¬ 
ables  were  used  (see  Figures  10  to  13).  The  r2  value  is  omitted  in  Figure  15. 
Since  there  is  no  y-intercept,  an  r2  value  would  represent  the  proportion  of  ex¬ 
plained  variability  about  the  origin.  This  value  cannot  be  compared  directly  to 
the  r2  value  computed  when  the  y-intercept  is  included. 


Predicted  vs.  Observed  Values  of  Understory  Stem  Density 


Predicted  Values  per  hectare 


Figure  15.  Predicted  vs.  Observed  Values  of  Understory  Stem  Density. 
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6  Summary 

Overstory  and  understory  tree  measurements  were  gathered  at  56  locations  in 
temperate,  mid-latitude  forest  sites  in  Central  Virginia.  Plots  were  selected  to 
provide  samples  from  the  full  spectrum  of  forest  structural  types.  Tree  height, 
average  crown  diameter  and  diameter  at  breast  height  (d.b.h.)  were  measured 
for  over  4300  overstory  and  understory  stems  on  10-m  radius  plots.  Other  plot 
information  included:  plot  slope  and  aspect,  drainage,  soils,  canopy  closure,  and 
types  of  vegetation  that  formed  the  herbaceous  layer.  This  information  was 
transferred  into  a  statistical  program,  and  simple  and  multiple  regression  were 
performed  to  find  the  best  overstory  variable(s)  to  use  in  predicting  understory 
stem  density. 

Overstory  mean  crown  diameter  proved  to  be  the  best  predictor  of  understory 
stem  density.  The  S.E.E.  for  the  simple  regression  of  overstory  mean  crown  di¬ 
ameter  vs.  understory  stem  density  was  520  stems  per  hectare.  This  translates 
into  95  percent  of  all  understory  stem  density  estimations  to  be  within  +/-  1040 
stems  per  hectare.  The  corresponding  coefficient  of  determination  (r2),  however, 
was  0.56.  This  indicates  that  only  slightly  above  one-half  of  the  variability  in 
understory  stem  density  could  be  explained  by  overstory  mean  crown  diameter. 
Multiple  regression  proved  to  be  a  less  accurate  predictor  (S.E.E.  =  573.5  stems 
per  hectare)  than  was  simple  regression. 
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7  Conclusions  and  Recommendations 

There  is  a  relationship  between  selected  overstory  parameters  and  understory 
stem  density.  When  graphed,  these  relationships  demonstrate  rather  structured 
patterns  that  are  logical  and  explainable.  The  standard  error  of  the  estimate 
(S.E.E.),  a  measure  of  the  dispersion  of  the  observed  values  about  the  regression 
line,  is  somewhat  excessive  and  the  coefficient  of  determination  (r2)  is  less  that 
desired.  This  may  be  due,  in  large  part,  to  the  fact  that  only  56  data  points  were 
available  for  analysis.  The  authors  were  unable  to  breakdown  this  data  set  fur¬ 
ther  by  species,  ranges  of  height,  crown,  and  d.b.h.  or  other  salient  characteris¬ 
tics.  Once  broken  into  smaller  subsets,  there  were  just  too  few  data  points  with 
which  to  derive  meaningful,  stable  equations.  This  was  especially  true  when 
multiple  regressions  were  performed  (e.g.,  different  sets  of  predictor  variables 
were  generated  when  forward,  backward,  and  standard  stepwise  regression  ap¬ 
proaches  were  used).  The  researchers  were  forced  to  use  the  data  set  as  a  whole, 
mixing  pure  stands  with  mixed  stands,  coniferous  with  deciduous,  etc. 

It  is  also  interesting  to  note  that  certain  variables  that  have  high  correlations 
with  vegetation  species  and  forest  density  at  other  geographic  sites  have  no  ap¬ 
parent  correlation  in  the  Virginia  piedmont.  For  example,  slope  and  aspect  have 
a  noticeable  relationship  to  species  and  density  in  most  mountainous  areas. 
However,  these  variables  correlated  quite  poorly  with  the  overstory  and  under¬ 
story  attributes  gathered  as  part  of  this  investigation. 

Additional  data  need  to  be  obtained  to  allow  analysis  by  species,  ranges  of  tree 
structural  measurements,  plot  characteristics,  etc.  More  work  needs  to  be  initi¬ 
ated  focusing  on  the  spatial  variability  of  the  understory.  Furthermore,  overstory 
and  understory  data  need  to  be  gathered  at  other  geographic  and  climatic  loca¬ 
tions  to  provide  an  indication  of  the  utility  of  any  developed  models  outside  of  the 
midlatitude  deciduous  forest. 

If  the  developed  models  are  eventually  deemed  successful,  then  research  needs  to 
be  conducted  on  the  overall  accuracy  of  the  image  analyst  in  measuring  the  re¬ 
quired  input  variables. 


40 


ERDC/TEC  TR-01  -1 


Bibliography 


Avery,  T.E.,  and  H.E.  Burkhart,  1994.  Forest  Measurements,  4th  ed.  New  York,  NY:  McGraw-Hill. 

DMA,  1982.  Defense  Mapping  Agency  Product  Specifications  for  the  Hard  Copy  Tactical  Terrain 
Analysis  Data  Base,  1”  ed.  Stock  Number:  SPECX050TTADB. 

HQDA,  1990.  Terrain  Analysis,  FM  5-33.  Washington,  DC:  Department  of  the  Army. 

Krause,  P.E,  H.B.  Puffenberger,  and  M.V  Campbell,  1999.  Structural  Relationships  of  Selected  Tree 
Species  at  Several  Mid-Latitude  Deciduous  Forest  Sites  in  Virginia,  Technical  Report  (TR)  TEC- 
0124.  Alexandria,  VA:  Topographic  Engineering  Center. 

Smith,  D.M.,  1962.  The  Practice  of  Silviculture.  New  York,  NY:  John  Wiley  &  Sons. 

SPSS,  1997.  TableCurve  2D  [Computer  Program].  Chicago,  IL:  SPSS,  Inc. 

StatSoft,  1999.  Statistica  for  Windows  [Computer  Program].  Tulsa  OK:  StatSoft,  Inc. 

Wenger,  K.F.  (Ed.),  1984.  Forestry  Handbook,  2nd  ed.  New  York,  NY:  John  Wiley  &  Sons. 


ERDC/TEC  TR-01-1 


41 


TEC  Distribution 


Chief  of  Engineers 

ATTN:  CEHEC-IM-LH 

Engineer  Research  and  Development  Center  (Libraries) 

ATTN:  ERDC,  Vicksburg,  MS 

ATTN:  Cold  Regions  Research,  Hanover,  NH 

ATTN:  Construction  Engineering  Research,  Champaign,  IL 

Fort  Huachuca,  AZ 

ATTN:  Information  Resource  Center 

U.S.  Army  Missile  Command 
ATTN:  AMSMI-RD-CS-R 

U.S.  Military  Academy,  West  Point,  NY 

ATTN:  Dept,  of  Geography  and  Env.  Engineering 


7 

3/01 


+12 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  mamtaniing  the 
data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing 
this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202- 
4302  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently 
valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS.  .  .  ..  - 

1 .  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

03-2001  Final 

3.  DATES  COVERED  (From  -  To) 

10-1999  to  07-2000 

4.  TITLE  AND  SUBTITLE 

Estimating  Understory  Stem  Density  from  Overstory  Structural  Characteristics 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

Paul  F.  Krause,  Michael  V.  Campbell,  and  Harry  B.  Puffenberger 

5d.  PROJECT  NUMBER 

ILIR 

5e.  TASK  NUMBER 

0060L5 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Topographic  Engineering  Center  (TEC) 

ATTN:  CETEC-TR-G 

7701  Telegraph  Road 

Alexandria,  VA  22315-3864 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

ERDC/TEC  TR-01-1 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Topographic  Engineering  Center  (TEC) 

7701  Telegraph  Road 

Alexandria,  VA  22315-3864 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

13.  SUPPLEMENTARY  NOTES 

Copies  are  available  from  the  National  Technical  Information  Service,  5285  Port  Royal  Road,  Springfield,  VA  22161. 

14.  ABSTRACT 

Current  airborne  imagery  provides  a  fairly  good  mechanism  for  the  measurement  of  such  forest  elements  as  tree  crown  diameter  and 
tree  liftiglit  The  nature  of  the  understory,  however,  can  only  be  subjectively  estimated  by  the  analyst.  Understory  density  and  species 
composition  can  have  an  impact  of  such  military  activities  as  cross  country  mobility,  cover  and  concealment,  line-of-sight,  and  bivouac 
locations.  This  work  involves  the  development  of  preliminary  models  to  estimate  understory  stem  densities  through  the  use  of 
measured  overstory  structural  characteristics  that  can  be  derived  by  an  image  analyst.  Over  4300  overstoiy  and  understory  stems  were 
measured  at  56  sites  at  seven  forested  locations  in  central  Virginia.  Measured  parameters  included  overstory  and  understory  tree 
height,  average  crown  diameter,  diameter  at  breast  height  (d.b.h.),  and  species.  Overstory  tree  locations  within  the  plot  (distance  and 
azimuth)  were  computed  and  plot  characteristics  were  gathered  (slope,  aspect,  ground  cover,  drainage,  soils,  etc.).  Preliminary 
equations  were  developed  using  simple  and  multiple  regression  to  estimate  understory  stem  density  from  these  selected  overstory  stem 
characteristics. 

15.  SUBJECT  TERMS 

tree  measurements  forest  understory  statistical  models 

forest  overstory  forestry  regression  analysis 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 

OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Paul  F.  Krause 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE 

Unclassified  Unclassified  Unclassified 

SAR 

42 

19b.  TELEPHONE  NUMBER  (in¬ 
clude  area  code) 

(703)  428-6802 

Standard  Form  298  (Rev.  8-98) 


Proscribed  by  ANSI  Std.  239.18 


Destroy  this  report  when  no  longer  needed.  Do  not  return  it  to  the  originator. 


